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We study some analytical properties of the solutions of the non perturbative renormalization group 
flow equations for a scalar field theory with Z2 symmetry in the ordered phase, i.e. at temperatures 
below the critical temperature. The study is made in the framework of the local potential approx- 
imation. We show that the required physical discontinuity of the magnetic susceptibility x(M) at 
M — ±Mo (Mo spontaneous magnetization) is reproduced only if the cut-off function which sepa- 
rates high and low energy modes satisfies to some restrictive explicit mathematical conditions; we 
stress that these conditions are not satisfied by a sharp cut-off in dimensions of space d < 4. 

By generalizing a method proposed earlier by Bonanno and Lacagnina ( Nucl. Phys. B 693 
(2004) 36.) to any kind of cut-off we propose to solve numerically the renormalization group flow 
equations for the threshold functions rather than for the local potential. It yields an algorithm 
sufficiently robust and precise to extract universal as well as non universal quantities from numerical 
experiments at any temperature, in particular at sub-critical temperatures in the ordered phase. 
Numerical results obtained for the ip 4 potential with three different cut-off functions are reported 
and compared. The data confirm our theoretical predictions concerning the analytical behavior of 
X (M) at M = ±M . 

Fixed point solutions of the adimensionned renormalization group flow equations are also obtained 
in the same vein, that is by solving the fixed points equations and the associated eigenvalue problem 
for the threshold functions rather than for the potential. We report high precision data for the odd 
and even spectra of critical exponents for different cut-offs obtained in this way. 
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I. INTRODUCTION 

During the last twenty years the non-perturbative approach to the renormalization group (RG) originated by 
Wilson fl], Q has been the subject of a revival in both statistical and quantum field theory. Two main formulations of 
the non perturbative renormalization group (NPRG) have been developed in parallel to study a system at equilibrium 
at, or near to, criticality ; for instance, in the simplest case, a system described, at microscopic scale A (in momentum 
scale), by an action 5a [<p] where if is a scalar field. In the first approach, one realizes a continuous RG transformation 
of the action <S/c[y] from k = A to k = and, a priori, no expansion with respect to whatsoever small parameter being 
required. At scale- k (in momentum space) the high energy modes tp q , \q\ > k, are integrated out in the "Wilsonian" 
action Sk[tp] which is a functional of the slow modes if q , \q\ < k. This "coarse-graining" operation requires the 
implementation of some cut-off of the propagator, either sharp or soft, aiming at separating slow (\q\ < k) and fast 
(\q\ > k) modes. This "coarse-graining" process is devised in such a way that all the Wilsonian actions Sk[<p] yield 
the same physics in the infra-red (IR) limit (q —> 0) . The flow of Sk [<p] , from the microscopic scale k = A to the 
macroscopic scale k = 0, is governed by the Wilson-Polchinski equation in case of a smooth cut-off and the 

Wegner-Houghton 0, Q equation in case of a sharp cut-off. Powerfull approximation schemes has been devised to 
obtain approximate, non-perturbative solutions of these equations; for a review see Ref. [f|. 

The second formulation, called the "effective average action" approach, was developed after the seminal works of 
Nicoll, Chang and Stanley for the sharp cut-off version [7], [|f and Wetterich, Ellwanger and Morris for the smooth 
cut-off version This method implements on the effective average action Tk[4>] - roughly speaking the Gibbs 

free energy of the fast modes <p q , q > k of the classical field <p = (ip) - rather than on the Wilsonian action Sk[<p] the 
ideas of integration of high-energy modes that underlies any RG approach. The flow of results in equations which 
can be solved under the same kind of non perturbative approximations than those used for the Wilson-Polchinski 
or Wegner-Houghton equations. The main advantage of this more recent formulation is that it gives access to the 
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RG flow of physical quantities, i.e. the Gibbs free energy Tk[4>] and its field derivatives, the vertex functions, rather 
than to such an abstract object as the Wilsonian action <Sfc[y]. Quite remarkably, the same kind of formalism was 
developed in a less elaborated language by Reatto et al. in the domain of the theory of classical liquids [l3l-[i"5l] , many 
years before these recent contributions to statistical field theory. The relations between these corpora of works is 
discussed in ref. . Recent reviews and lectures devoted to Wetterich's approach are available [l?], [H| and should be 
consulted for a thorough discussion. Wetterich's version of the RG is in fact equivalent to that of Wilson- Polchinski 
as discussed in Ref. [12, Eil ■ 

In this paper we will adhere to Wetterich point of view and focus on the study of approximate solutions of the 
NPRG for a scalar field theory with Zi symmetry, at a temperature T below the critical temperature T c , i.e. in the 
"ordered phase" or "two phase region". For T < T c the system exhibits a spontaneous magnetization M =< ip > 
which can take any value in the interval (— Mq(T), M (T)) with an equal probability in the thermodynamic limit. 
Henceforth we shall reserve the name spontaneous magnetization for Mq(T). As a result the susceptibility x is infinite, 
or equivalently, the second derivative of the Gibbs free energy with respect to the magnetization 

r"(M) = x^ 1 = in the ordered phase. The potential T(M), strictly convex for \M\ > M (T" (M) = x~ X > 0) is 
affine for -M < M < M . 

The simplest non-perturbative version of the NPRG, i.e. the local potential approximation (LPA), yields indeed a 
convex free energy with a plateau as noticed in Refs fl3T4l5l UtI [20j . However it has been discovered by Reatto et al. 
that the analytic behavior of T (M) in the vicinity of ±Mo depends strongly on the choice of cut-off and dimension 
d. For a sharp cut-off (and in d = 3), T (M) is a continuous function of M, notably at ±Mo, and thus T (Mo+) = 
[l3T |. This is a serious flaw of the theory since obviously T (M) should be discontinuous at ±Mq, i.e. vanishing 
identically in the two phase region with a jump to a finite positive value outside the two-phase region corresponding 
to a finite susceptibility. In other words one should have T (Mo—) = and T (Mq+) > 0. However, with another 
choice of cut-off function proposed by Litim [21], it was shown by Parolla et al. in Ref. [l5[ that, at least in d — 3, 
r (M ) exhibits the correct discontinuity at M = ±Mq. In this work we extend these results to any kind of cut-off 
and dimensions d of space and obtain the mathematical properties to be satisfied by the cut-off function to obtain 
the required discontinuity of the susceptibility at M = ±Mq. The crux of the whole matter is that the RG flow of 
Gibbs free energy stops in the ordered phase at some finite RG time and that the solutions of the RG flow can thus 
be obtained as asymptotic stationnary solutions for well chosen dimensionned functions the behavior of which gives 
insight on the region M ~ Mq. 

These analytical results are then checked by numerical experiments in the case of a <p A potential where the RG 
flow equations are solved with a new algorithm which generalizes that devised by Bonanno and Lacagnina (20j to 
any kind of cut-off. The idea is to solve the RG flow for the threshold function rather than for the potential itself, 
or one of its derivatives. The partial differential equation for the threshold function is of quasi-linear parabolic type 
(rather than of non-linear parabolic type for the potential yielding huge numerical instabilities in the ordered phase) 
and can thus be solved with an arbitrary high accuracy. The approach to convexity of the Gibbs free energy is then 
achieved exactly with V (M) vanishing identically in the ordered phase, down to the smallest real available on your 
computer if wanted. Solving RG flow equations for the threshold function above T c is also possible of course, but 
of less interest. A good numerical precision can therefore be achieved making possible to extract precise universal 
as well as non-universal quantities from numerical experiments. Thorough numerical investigations have thus been 
undertaken with three different cut-off functions, i.e. the sharp cut-off, Litim's regulator and an exponential smooth 
cut-off, aiming inter alia at testing the theoretical predictions on the behavior of x(M) at M = ±Mo; a thorough 
comparison and discussion of the data provided by the different cut-off is also presented. 

This study is finally completed by solving the adimensionned flow equations asymptotically in the same vein as 
that used to solve the dimensionned equations, i.e. by solving the fixed point equations and the associated eigenvalue 
problem for the threshold functions rather than for the potential. Fixed point functions and critical exponents are 
then obtained for the three cut-offs with a very high numerical precision, notably for the exponential smooth cut-off 
non considered up to now; in the case of the sharp and Litim's cut-off we recover the results obtained previously 

The paper is organized as follows; in section |TT] we summarize briefly Wetterich formalism and the LPA approx- 
imation. The properties of the RG in the two-phase region are then explored in depth yielding the mathematical 
conditions to be fulfilled by the cut-off function in order to obtain the correct discontinuous behavior of T (M) at 
M = ±Mq. In section InTl we devise a new numerical algorithm in which a correct account of both initial and boundary 
conditions is given. This material is used to perform extensive numerical experiments aimed at testing the theoretical 
predictions of section |H] and at computing some universal and non- universal quantities of the (p 4 potential. Section HVl 
is devoted to a somehow new presentation and numerical study of the adimensionned flow equations. Estimations 
of the critical exponents of the <p A (or Ising) 3d model in the LPA and for three different cut-off are reported. In 
the case of Litim's and sharp cut-off our results are in perfect agreement with those of ref. [22|, [23| and in the case 
of the exponential smooth cut-off we provide the first high precision data for the even and odd spectrum of critical 
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exponents. We conclude in Section W\ 

II. NPRG FLOW EQUATIONS IN THE LOCAL POTENTIAL APPROXIMATION 

A. The model and the NPRG 

We consider a scalar field theory defined at scale A by it's microscopic action 

SaM = J x |(v^) 2 + ^a(^)| , (1) 

where cp x is a scalar field defined at some point x of the domain f2 which encloses the system, J = d d x (d = 
dimensions of space.), Z2 symmetry is assumed (i.e. V\ (<p x ) — V\(—(p x )) and we stick to the convention that 
brackets [. . .] denote functionals while parenthesis (. . .) are for functions. Below, for numerical applications, we shall 
restrict ourselves to the (p 4 potential , i.e. a coarse-grained version of the Ising model with V\ ((fx) = ^\ r f x + jiff^x 
( g > 0), and will denote by r c the value of parameter r at the critical point in the LPA approximation. Finally it is 
understood that A acts as an ultra-violet (UV) cut-off in integrals defined in Fourier space, i.e. j q = d d q/(2ir) d . 

Following Wetterich [l7| we introduce a family of related fc-models depending on an index k in momentum space 
with < k < A. The bare action Sk [<p] of the k— model - not to be confused with the Wilsonian action briefly alluded 
to in the introduction- is defined by 

Sk M = S A [ip] + \p ■ !Z k ■ <p , (2) 

where ip ■ TZk • P> = J p q <Pp T^kip, l) is a massive term aimed at separating high and low energy field modes. TZkip, q) 
is called the cut-off or regulator function. Translational invariance implies lZk{p, q) — Rk(.q 2 )(^) d 5 d {p + q) and the 
general behavior of the cut-off function is Rk{q 2 ) — Zkk 2 (l — e (<?, k)) where Zk is some positive prefactor accounting 
for field-renormalization and Oe(q, k)) is some smooth approximation of the Heaviside step function Q(q — k) so that 
when e — > then & e (q, k) — > <d(q — k). The cut-off functions must have the following generic behavior 

(i) when k = 0, Rk=o(q 2 ) = identically (Vg) and the original model is recovered. 

(ii) when k = A, Rk=A(l 2 ) — Z^K 2 is sufficiently "large " so that all the fluctuations are frozen by the massive 
term, i.e. the mean-field approximation becomes nearly exact. 

(hi) when < k < A, Rk(q 2 ) is a decreasing function of q which tends rapidly to for \q\ > k, thus the rapid modes 
are unaffected by the massive term. On the contrary the slow modes have a large mass which decouples them 
from the fast modes. 

The physics of the k-systems is encoded in their partition functions 

Z k [J] = J V^cxp (-Sk [<p]+J- <p) , (3) 

where the functional measure Dip integrates over all the field modes cp g (\q\ < A) and J ■ (p = f J x ip x (J source). 

At scale "k" the slow modes <p q , \q\ < k, are frozen by their mass Rk(q 2 ) ~ Zkk 2 and Wk [J] = log Zk [J] can therefore 
be interpreted as the generator of the connected Green's function with an IR cut-off. It is a convex function of the 
source J and it's Legendre transform defined as 

r fc [0] =sup(J-cj > -W k [J}) , (4) 
J 

is also a convex functional of the order parameter (f> —< ip >. r& [<fi] is the " true" Gibbs free energy of the k-system. 
However it proves convenient to introduce and consider rather the effective average action of Wetterich 

r fc [0] = T k [ct>]-^<p-ii k -(t>. (5) 

By contrast with T k [4 1 ] this functional of the field 6 is non-convex but has simple limits which follow from the 
mathematical properties of the cut-off function TZk [17| . i-e. : 
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• when k = A no fluctuations has been integrated out. and we will suppose that Tk=\ [4>] = $A [4>] (mean field 
approximation) . 

• when k = all fluctuations have been integrated out and Tk=o [(f)] = T [(f)] is the Gibbs free energy of the model 
under study. 

Therefore when k decreases from its initial value A to k = all the fluctuations of the field are progressively 
integrated out and the effective action [(f)] flows from the bare action 5a [(f)] (i-e . the mean-field approximation for 
r [(f)]) to the Gibbs free energy T [(f)] of the model. The flow equation reads as [13] 

d k T k [4>] = \J d k TZ k {-q, q) {K k + f { k 2) }-\-q, q) , (6) 
where the inverse in the r.h.s. of Eq. © has to be understood in the operator sense, T k 2 \x : y) = S 2 Tk[(f)]/S(f) x S(f)y 

— (2) (2) 

is the vertex function of order 2 and T y k '(p, q) denotes is 2d dimensional Fourier transform. Note that T k ' depends 
functionally on the field (f> and the equation for Tk [(f)] is thus not closed. One has to resort to approximations to solve 
the flow. 



B. The local potential approximation 

A simple, non trivial way to tackle with Eq. ([6]) is to restrict the functional space to functionals of the form 

(LPAansatz) T k [0] = J |i(V^) 2 + C/ fc (^)| , (7) 

which constitutes the popular local potential approximation (LPA). Note that in this scheme there is no ficld- 
renormalization so that Zk = 1. We also stress that the local potential Uk((f) x ) is a function of the local field 
<j) x at point x, not a functional. Combining Eqs. ([6]) and ([7]) one then obtains a closed partial differential equation 
(PDE) for the Uk which reads as 

where the potential Uk(M) is evaluated for a uniform magnetization <fi x = M and U k (M) = d 2 Uk(M)/dM 2 denotes 
its partial derivative with respect to M. The PDE Eq. (jHJ is non-linear but however quite easy to solve numerically 
above the critical point (r > r c , in the case of the ip potential). Difficulties arise in the ordered phase. To understand 
why, we rewrite Eq. (J5J under a more convenient form. We define the dimensionless variable t — — log(fc/A), the 
RG time chosen to increases from to +oo when k decreases from A to 0, the dimensionless variable y = q 2 /k 2 
and the dimensionless functions of "y" : r(y) = Rk(q 2 )/q 2 (so that Rk(q 2 ) = yr(y)k 2 ), s(y) = —2y 2 r'(y), and 
t{y) = y{l + r{y)) so that 

d t U k (M) = -2v d k d C(U' k \M)/k 2 ) , (9) 
where v^ 1 = 2 d+1 n d r(d/2) and the threshold function C (according to Wetterich terminology) is defined as 

£: ' r »«-''(»)*. do) 

2 Jo t(y)+u 

This function C(ui) plays a central role in the mathematical analysis of Sec. Ill CI and we need first to study its 
relevant analytical properties. For reasonable cut-off functions the two functions s(y), and t(y) are positive for y > 
so that C{uj) is a decreasing function of it's argument. In any reasonable case the function is defined on the interval 
]l<jo,+oo[ where is the largest pole of the integrand in the r.h.s. of Eq. (|10p (for any reasonable choices of the 
cut-off r(y) there is in general only a single pole). Therefore C(uS) decreases from +oo to when u> increases from ujq 
to +oo. Finally note that U k (M)/k 2 is indeed a dimensionless quantity which gives sense to Eq. ([9]). 

Before we tackle the general case let us consider two important cases. First, the (ultra) sharp cut-off defined by the 
regulator Rk(g 2 ) — (3k 2 ( 1 — 9 (q 2 — k 2 )) where the parameter j3 — > +oo. In that case one shows that £(ui) =— ln(l+a;) 
so that the flow equation reads as 

d t U k (M) = 2v d k d ln(l + U' k \M)/k 2 ) . (11) 
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A second important case that we shall consider is Litim's regulator Rk{q 2 ) — (k 2 — q 2 )(l — 9(q 2 — k 2 )) which yields 
C(u) = (2/d)/(l + uj) and the flow equation [H| 

8 > u *w = -^TTumw- ■ (12) 



TABLE I: Constants characterizing the exponential smooth cut-off threshold function in the case a = 6. 



2/o = 


2.36806229624554475822396946574 


wo = 


-3.83637249545350290077137812782 


A'«, = 


34.2450749455029367943794121780 


K l = 


26.7493978693695013104149308428 


K 2 = 


-123.598329790226947089342622508 


K 3 = 


600.184183291720248551501270440 


Ci = 


26.7493978693695013104149308428 


C 2 = 


-4.62060231762294359247562464639 


C 3 = 


.0406488093691634129962437962261 



Let us consider now a regular smooth cut-off function. We have retained the exponential regulator widely used in 
recent numerical studies (see e. g. ref 

R k (q 2 ) = aq 2 /(e X p(q 2 /k 2 ) - 1) (13) 
that is, written in reduced form r(y) = a/(exp(y) — 1), where a is some positive parameter. It is easy to see that 

(i) When a > 2 the function t(y) = y(r(y) + 1), defined for y £ [0, +oo[, exhibits a single minimum at some yo > 
and t'(y ) =0. 

(ii) When a = 2 the minimum of t(y) occurs precisely at yo = and t'(0) = 

(hi) When < a < 2 the minimum of t(y) on the interval (0, oo) is located at yo = but t'(0) > 0. 

The choice a = 6 yields the better critical exponents according to the authors of reference [24| so we will retain this 
value; we are thus in the case (i) where the minimum of function t(y) is located at some non-zero value yo- We 
introduce ujo = — min o<y<oo t(y) = —t(yo). The threshold function uj t— > C(u>) is thus a monotonously decreasing 
function defined on the interval ]u>o, oo[. Let us precise now its asymptotic behavior at each boundary of the interval. 

Wetterich et al. [17j have shown how to obtain the behavior of C[ui) for w — > cjo+; one expands t(y) around its 
minimum, i. e. t(y) = —uia + ^2^y 2 + 0{5y i ), where 5y = y — yo and recognize the fact that, at the leading order, 

2 J-oo uj-ujo + t 2 dy ! 
~ K (lo - uo)- 1/2 , (14) 

Ko = 2vr 2 ' (15) 

and t 2 = t"(y )/2. 

Note that in Litim's case one has loo = — 1, t(y) = 1; therefore the previous analysis breaks down and C(u>) diverges 
as ~ (oj — wo) 1 rather than as ~ (u> — cj ) -1 ^ 2 - In the case of the ultra sharp cut-off u>o = —1 but the function t(y) 
is not well defined yielding a logarithmic divergence of C{uS) as oj — > loq = — 1. 

On the other hand an asymptotic behavior of C(u>) for uj — > oo is readily obtained from (|10[) and reads 

r( \ K 2 K 3 

C(oj) = — + — + — + ... 
uj cu z u; J 

K n = ( —!— / dy y^-h^tiyr- 1 . (16) 



where 




Since uj > L = C(uj) is bijective from ]wo,oo[ to ]oo,0[ it can be inverted and we denote by 
This function is defined on the interval ]0, oo[ where it decreases from +00 to ojq. 
For L — > we infers from (|TB)l that 



j(L) « ^ + C 2 + C 3 L + . . . , (17) 



L 



where 



and for L — > 00 



K2 n _ K 3 - K% /K\ 



C 1 =K 1 ; C 2 = ■ C 3 = T J< , (18) 



K 2 

oj(L)=w + j± + ... . (19) 

The values of the constants yo 7 uJo, Kq, K±, K2, if 3, Ci, C2, C3 are determined numerically; they are resumed in Table U 
in the case a = 6. For numerical applications the functions uj{L) or £(w), as well as their derivatives if required, have 
been fitted by polynomial expressions taking into account their asymptotic behaviors. 



C. Flow equations for the threshold functions 

As suggested first by Bonanno and Lacagnina [20| (however in a more restricted context) it is convenient to perform 
the change of variables (M, k, U) (M, k, L) with 

L k (M)=£(U>,!(M)/k 2 ) . (20) 

The one to one mapping uj ■(=> C insures the mathematical equivalence of solving the RG flow equations either fo r U'J 
or for ij.. However, for numerical reasons the RG flow equation for L k is much easier to solve than that for Uk |2fj| . 
The RG flow equation for Lk is readily deduced from that for Uk {cf. ©) : 

L'l = --—uj\L k ) -l± + —uj(L k ) . (21) 
2vd at Vd 
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FIG. 2: Litim's approximation : universal function L},{M)/Lk(M = 0) as a function of x = M 2 /M(0) 2 for different value of 
the scale k m in < k < A. Universal curve (I26p : red line ( k = fcmin = 10 -6 ). The </? 4 potential is used at scale A = 10 with 
g = 1.2 and r = -0.35, M max = 8. 



Eq. (j2Tj) is a quasi-linear parabolic PDE which can be studied analytically in some depth, notably in the two phase 
region, and for which, in addition, the mathematicians have provided us with robust and efficient solvers in view of a 
numerical study. For the exponential smooth cut-off (|13l) the function w(Lfc) and its derivatives which enter Eq. (|21[) 
must be determined numerically. It is worthy to write down explicitely these flow equations 

• for the sharp cut-off regulator (C(uj) = — m(l + to) and oj(L) = cxp(— L) — 1) 

L'l = — exp(-L fe )— * + (exp(-L fc ) - 1) . (22) 

2v d at v d 

• and for Litim's regulator (£(w) = (2/d)(l + lu) and u>(L) = 2/ '(dLk) — 1) 

+ (23) 
v d dL% dt v d dL k 



D. Behavior of the flow in the ordered phase 

Below T c , a spontaneous magnetization Mo(r) can settle in the system in the absence of an external magnetic field 
J. More precisely, in the thermodynamic limit, any magnetization — Mo(r) < M < Mo(r) is likely to settle with the 
same probability. Roughly speaking we have a coexistence region of two (or more) magnetized phases ±Mo analogous 
to a liquid-vapor coexistence. Therefore T"(M) = dJ/dM oc r^ 2 ^(0) = from which it follows that, in the limit k — >• 
the denominator of the r.h.s. of Eq. (J8)) tends to zero. Therefore the threshold function Lk(M) should become large 
and positive as k — » and finally it should diverge to +oo at k — for any magnetization in the two phase region 
-M < M < M Q . 

With the hypothesis that Lk is large and positive, ui(Lk) ~ ujq as follows from (|T^)) ; then the flow equation (|2"Tj) for 
Lk simplifies to 

1,2-d 

L'l = Wo , (24) 

Vd 




FIG. 3: Litim's approximation : from right to left and bottom to top : local potential Uk(M), the derivatives U k (M) (magnetic 

field) , U k (M) (inverse susceptibility), and U k (M) for different values of k in the range fe m i n < k < A, with fc m in = 10 . The 
ip 4 potential is used at scale A = 10 with g = 1.2, and r — —0.35 (M max = 8 ). Red lines : result at k = fc m i n . 



which can be integrated to give 

L fc (M) = uj (Mo(fc) 2 - M 2 ) . (25) 

for — Mo(fc) < M < Mo(k) where Mo(k) is the constant of integration (depending on k) of Eq. (IM1) . Clearly Mo(fc) 
can be interpreted as a precursor of the magnetization of the system at scale k. Note that M (k) can be obtained 
numerically from the numerical solution of PDE (|21[) from the value of L k (M = 0) = — ujoMo(k) 2 k 2 ~ d / (2vd) at the 
origin. As expected one finds that for d > 2 the function L k (M) diverges to +oo (recall that uj n < 0) when k —> as 
oc k 2 ~ d for any magnetization M in the two phase region. 

Outside the coexistence region -i.e. for \M\ > Mo(fc)- one expects a finite compressibility and thus L k (M) tends 
to as ~ K\k 2 /U k (M) with k, as follows from (fT6]l . This behavior of L k (M) is exemplified in figure [T] (details 
concerning the numerical procedure to solve (f2Tj) will be given in Sec. III). 

One of the consequences of Eq. (|2"5)) is the universal behavior : 

L k (M)/L k (0) = l-x with x = (M/M (k)) 2 for | x \< 1 

L k (M)/L k {0) = for | z |> 1 . (26) 

The function L k (M) / L k (0) of the argument x = (M/Mo(k)) 2 is thus universal as k — > , i.e. it is independent of 
the thermodynamic state (provided the temperature satisfies r < r c ), the regulator R k (q 2 ), the value of ujq, and the 
dimension of space d > 2. This remarkable features are illustrated in figure f2J 

One can of course infer the behavior of the potential U k {M) and its derivatives from that of L k and one has, for 
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-Mo(Jfc) < M < M (Jfe) 

£/£(M) = fc 2 w 

17* (M) = fc 2 w M 

l/ fe (M) = - k 2 uj M 2 + constant , (27) 

where we have noted that U'AM = 0) = for a Z2 symmetry; see figure [3] for an illustration. In the special case of 
the sharp cut-off and Litim regulator recall that one has to set ujq = — 1 in the above equations. Another interesting 
consequence of this simple behavior for the potential and its derivatives is the expression of the "true" susceptibility 
of the k-system (cf Eq. (|5j)) which reads as 

-^r fc (M) = u" k {M) 



dM 



= U k (M) + yr(y)\ y ^ 

= k 2 (t(0)- min t(y)) , (28) 

and is thus positive at any k ^ 0. However the strict convexity of Tk(M) is preserved at each step of the RG flow in 
the LPA approximation only if a > 2 as follows from the discussion after Eq. (fT2"j). Note that in that case Eq. I|28p 
d 2 - 

becomes ^Tfc(M) = (a + u>o)k 2 . In the case a < 2, strict convexity could be violated and the validity of the LPA 

dM z 

is no more guaranteed. 

The simple behavior described by Eqs. ([2~T]) is independent of the regulator and has been known from long [LH [20|. 
However Parolla et ah (l3l - fIH were the first to remark, in a more refined discussion, that the behavior of the 
susceptibility for M — > ±Mo(fc) depends strongly on the asymptotic behavior of C(ui) for oj —> ujq+. For the sharp 
cut-off regulator they found that U k (M) is continuous and equal to zero at M = Mq (for d < 4) (see Ref. [l3j]) while 
for the Litim regulator they found that, for d > 2, TJ V (M) exhibits a discontinuity from a finite positive value at 
M = Mq+ to a zero value at M = Mo— ( cf Ref. This discontinuity of the susceptibility \ jumping from an 

infinite value in the two phase region to a finite one outside, is of course the expected physical behavior for \. In 
other words, in the case of the sharp cut-off regulator, the spinodal and the binodal curves merge in a single curve, 
which is a serious flaw in the theory (see figure 2]) . 

We now extend the analysis of Parolla et ah to the quite general case : 

C(u) = K ° , (29) 

where v > is an arbitrary exponent, not to be confused with the critical exponent of the correlation length! This 
expression includes notably the case of the exponential smooth cut-off considered in this paper - with v = 1/2, cf 
Eq. (fl4|) - since only the singularity of C(u>) as uj — > loq+ is relevant to study the behavior of U k (M) as M -> ±Mq. 
As follows from (|25|) . when k —> L k (M), diverges therefore to +00 for \M\ < Mo for d > 2 and tends uniformly to 
outside the interval ] — Mo, Mq[ (see fig. [1]). This behavior is too "drastic" to discriminate the analytical properties 
of U (M) in the vicinity of ±M . Borrowing an idea of Rcatto et ah we introduce the following new function of the 
field 

F(k,M) = - — — - — - = K^k- 2v L k {M) . (30) 

In the two phase region U k (M) — » as k — > and thus F — > +00, while outside < U k=0 (M) < 00 yielding F to 
reach a finite value F at k = 0. It follows from ([25]) and ([30]) that for k -> and \M\ < M (k) 

F(k, M) = k 2 - d -' 2u uoo{x 2 -xo{k) 2 ) 
M 



Mo(fc) 

zo(fe) = = ■ (31) 

V2fdA 

We are led to introduce the new variable 

z = (x- x (k))k 2 - d - 2u (32) 
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FIG. 4: Uu. (M) with sharp (green solid circles), Litim's (black solid circles) and exponential smooth cut-off (red solid circles) 
in d — 3 dimensions of space. In all cases : ip 4 potential at A = 10 with g — 1.2, M max = 8, SM — 2.10 -6 , and fc m i n = 10 -8 . 
Litim's cut-off : r = -0.197032, sharp cut-off : r = -0.297286, exponential smooth cut-off : r = -0.412195. 



such that, at fixed z and for k — > 0, one has 

F(k,M) tt2u x (k)z , (33) 

valid if 2v > — d + 2, i. e. v > if d > 2. As quoted pleasingly by Parola et al.[lj£, UU this variable z allows us to 
"zoom" in the region x ~ x (k). Inside the two-phase region z — > — oo and we have the asymptotic behavior (|33[) of 
F, while for z — > oo, i. e. outside the the two-phase region, the function F should reach a finite value Fq. To avoid a 
proliferation of notations we still write F{k, z) the function F(k, M) expressed in terms of its new variables (fc, z). 

The dependence of the reduced spontaneous magnetization xo(fc) on k as k — ¥ will prove of great importance in 
our analysis. A priori it should be reasonable to assume a k-dependence as 

M (fc) - M - ak 2 with a > (34) 

for sufficiently small k's, since Mg(k) should be an analytical function of vector k. Parameter a > describes 
the displacement of the precursor of the spontaneous magnetization Mo(fc) = \/2vdKoXo{k) with scale k. With this 
reasonable assumption we can finally wrote the RG flow equation : 



z = kf- d -^{x-x + ak 2 ) 
k*- d - 2v -^F(k,z) = 2fcW ' ' 



dz 2 y ' ' u v F 1 /^+ 1 {k.,z)) 

x [kd k F(k, z) + (2-d- 2v)zd z F(k, z) + 2ak i ~ d ' 2u d z F{k 1 z)] . (35) 

We now analyse the asymptotic behavior of the stationary solutions (kdkF(k, z) — 0) of EDP ([33]) at z — > ±oo in 
the limit k — > 0. This is justified since below r c the flow stops at some finite k. These solutions are similar to scaling 
solutions at a fixed point; however F(k, z) and z have dimensions and we cannot thus speak of scaling solutions stricto 
sensu. Obviously, in order to find these solutions, which we cannot refrain to christen as fixed point solutions despite 
our previous remarks, we have to discriminate our study according the sign of the exponent 4 — d — 2v in Eq. (|35[) . 



0.0006 



£ 0.0003 




0.0478 0.048 

M 



FIG. 5: t/ fcm . n (M) for the simplified smooth cut-off (i. e. v = 1/2 in Eq. p9)ll in the vicinity of M . </? model with A = 10, 
g = 1.2, fc min = 10~ 8 , M Max = 8, d = 3.25 and r = -0.1545. Color code for the open circles (if present) : black : SM = 2.10" 3 , 
red : 5M = 2.10" 4 , green : SM = 2.10~ 5 , blue : <5M = 2.10~ 6 . All the points of the grid of field values are reported, dashed 
lines are only guide-lines for the eyes. 



4 - d - 2v < 



Neglecting terms which tend to when k — > in Eq. (|35j) one finds the implicit fixed-point equation 

d 



dz 2 y 1 v F 1 /^ 1 (z)dz 



(36) 



where F(z) = F{k = 0, z). This is an autonomous differential equation which can be solved analytically : 

-cz + d = F Q F(F/F ) , (37) 

where c > 0,d are integration constants, Fq = (2ajc) u and T{X) a primitive of f(X) = 1/(1 — X^ 1 ^). We define no 
as the integer part of v, 5 as its decimal part, i. e. v = uq + 5, and we have 



F(X) = X + - {X 5 '» - 1) + v (1 - 5 n0)1 ) ^ 



v — n 



vX^-^l^iX- 1 !",!,!- 5) , 



(38) 



where the third term in the r.h.s. of the equation is present only if v > 2 and $(z,s,a) is the Lerch transcendent 
[25|,[26j]. This curiosity is defined on the disk \z\ < 1 (although it can be extended by analytical continuation to the 
cut plane C — [1, oo[). It is defined as 



(a + n) s ' 



0^0,-1,..., \z\ < l;\z\ = l,3?(s) > 1 



(39) 
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When v — no is an integer Eq. (|38[) reduces to the simple form : 



+ n ln(l -X- 1 /" ) . (40) 

This result was first obtain in Ref.[l5| in the case no = 1 (Litim case). From the known behavior of Lerch function, 
in particular $(z~ 1 / l/ , 1, 1 — S) ~ — ln(l — z^ 1 ^) for z — >• 1+ one deduces that 

- X for X -> +oo 

J'(X) - v\n(X 1/v - 1) for X -> 1+ (41) 

We are now in position to discuss the property of the fixed point solution Eq. (|37[) . In the two phase region, z — > — oo, 
the function F ~ 2ujqXqz diverges to +oo, thus F(F/Fq) ~ F/Fq, and one has asymptotically — cz ~ 2u>qXqz which 
fixes the value of the integration constant c = — 2wqXq, from which we deduces Fq = (— a / (wqXo)) 1 ' . In the one phase 
sector, z — > +oo we find that F — > Fq— reaches a finite value as fc — > 0. Returning to our old variables (k,M) we 
have shown that the fixed point solution satisfies : 

lim lim F(k, M) = < F Q < oo , 

M-s-M + fc-s-0 

lim lim F(k, M) = +oo , (42) 

M->Mo-M 

from which follows that for 4 — d — 2v < (and <i > 2) the susceptibility \ = M) v is discontinuous at M = ±Mq. 
See Fig.© for an illustration. 

2. 4~d-2v = 

We consider here the marginal case 4 — d — 2v = 0. It is an important one since it corresponds to the general 
smooth cut-off (v = 1/2, cf Eq. (jT4j) ) in d = 3. Neglecting terms which tend to when k — ► in Eq. (|35|) one finds 
the implicit fixed-point equation : 

^'W + f-^PT^wS'W- (43) 

This equation does not seem to have an analytical solution but, asymptotically : 

F(z) <~ 2uiqX()Z for z — > — oo , (44a) 
F(z) ~ F < +oo for z -> +oo . (44b) 

Equation (I44al) is required by the correct matching with the solution (1331) in the two-phase region while (|44bl) guar- 
antees finite susceptibility outside the magnetization curve. We conclude again in favor of a discontinuity of the 
susceptibility at M = Mq. See Fig.© for an illustration. 

3. 4-d-2u>0 

In this case the term in "a" which accounts for the displacement of the spontaneous magnetization Mo(fc) with 
scale k (cf Eq. can be discarded from Eq. (|3"5)) . The function F(k, z) is no more adapted to our discussion and 

we are led to a new change of variables to eliminate the relevant dependence on k from the RG flow equation. Let us 
introduce 

G{k,z) = k~' i F(k,z) , 
d-±-2v 

fx = — — v >0. 45 

v + 1 

By use of these variables the RG flow equation for G admits in the limit k — > a stationnary solution G(z) which 
satisfies the following differential equation : 

£ G(z) + I g^(z) { ^ g{z) + (2 " d 2 ^ Z §; G ^ = ■ ( 46 ) 
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FIG. 6: Same as in fig. [5] except that d — 3. and r = —0.1495. 



This equation cannot be solved analytically but it satisfies the boundary conditions 



G(z) oc z for z — > — oo , 
G{z) ~ z /»/(««-a-ai/) for 



(47a) 
(47b) 



Equation (|47a[) is required by the correct matching with the solution (1551) in the two-phase region, while (I47b|) implies 
a non-physical divergence of the susceptibility on the magnetization curve. Indeed it follows from (|47b[) and (|45l) that 
[/£ =0 = i?- 1 /" oc (M - M )™ (M > M since z -> +oo)) where the exponent 



is strictly positive if d > 2. Not unexpectedly, one recovers, in the limit v — > 0, the exponent w = (4 — d)/(d — 2) 
obtained by Parolla a/, in the case of the sharp cut-off [ljj. For the cut-off considered in this section the LPA 
seems to detect the spinodal rather than the coexistence curve and the exponent w can be interpreted as the exponent 
of the inverse compressibility x _1 (M) on the spinodal. Note that, however, it depends strongly on the non physical 
exponent v which characterizes the divergence of the threshold function as uj — > ujq. 

We conclude that for 4 — d — 2v > (and v > and d > 2) the LPA has the unphysical feature to give an infinite 
susceptibility at M = ±Mo- The isotherm U (M) is continuous with M and vanishes exactly in the two- phase region. 
See Figs. ([7]) and © for an illustration. 

The conclusion of this section is that, if the singularity of the threshold function u at lj = ujq+ is characterized by 
the exponent v (cf Eq. (1291) ) and if d > 2, then 

4— d 

• For < v < — - — the inverse compressibility U"(M) is continuous at Mo, which constitutes a severe flaw of 



4-d + 2v 



(48) 



zu = 



(v + l)(d - 2 + 2v) 



the theory. 



A-d 



< v the inverse compressibility U"(M) is discontinuous at Mq 



• For 
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FIG. 7: Same as in fig. [5] except that d = 2.75 and r = -0.1495. 

III. NUMERICAL INTEGRATION OF THE DIMENSIONNED RG FLOW EQUATIONS 

A. Algorithm 

In order to solve the EDP (j2lT) we must specify 

• (i) an arbitrary initial condition at t = (k = A) 

L A (M) = Li[M) VM e (-M max ,M max ) , (49) 
where M max is some maximum magnetization. 

• (ii) boundary conditions 

L fc (±M max ) = L b (k) Vfce(fc min ,A). (50) 

where fc m in is some minimum value of the scale k. 

• (iii) The functions L^M) and a priori arbitrary must be such that [30l | 

Li(±M max ) = L 6 (A). (51) 

Under these 3 conditions the RG flow equation for have a unique solution. The choice of Li and relies on 
physical grounds. At k = A the cut-off function R oc oo should diverge [13, EH such that Ua = V\. Thus, clearly, 
Li(M) = C(V A (M)/A 2 ) must be imposed as an initial condition even if A is not large enough so that Ra is not 
actually infinite in the mathematical sense. 

We propose to choose as boundary conditions the one- loop approximation for Lfc(±M max ). Most authors usually 
adopt free boundary conditions. At a large M max fluctuations are frozen and the one-loop approximation should be 
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FIG. 8: Same as in fig. [5] except that d = 2.5 and r = —0.157. 



pertinent. It is easy to see that at the one-loop level one has, up to an additional constant, independent of the field 



Ul Aoov (M) = 



^1-loop 



(M) 



V k (M) 

1 f, r q 2 + V^(M) + R k (q 2 ). 



from which it follows that 



5 fc ^- loop (Af) 



1 



A 2 



d k R k (g 2 ) 



(52) 



(53) 



(54) 



2 h q 2 + V^(M) + R k (qi) 
= 2v d k d - 1 C(V k {M)/k 2 ) 

Therefore the one- loop approximation for L k (M) reads as 

To summarize we impose 

• (i) initial condition : Li(M) = £(V^'(M)/A 2 ) for all - M max < M < M max . 

• (ii) boundary conditions Lb(k) = C(V k (±M max )/fc 2 ) for all < k < A. 

and we stress that the compatibility condition (jSTj) is obviously satisfied hence the RG flow equation can be solved 
safely at least from a mathematical point of view. 
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t 



FIG. 9: U k (M = 0)/(fc 2 |u;o|) as a function of the RG time "t" in the Sharp cut-off approximation. r c is obtained by dichotomy. 
Dashed lines : r > r c . Solid lines : r < r c . 

B. Numerical experiments 

The RG flow equation (|2Tj) , completed with the initial and boundary conditions of Section IIII A[ is a standard 
well-conditionned quasi-linear parabolic EDP which can easily be solved numerically. We adopted the fully implicit 
predictor-corrector integration scheme of Douglas and Jones j|32j which proved to be very efficient and precise. Douglas 
and Jones have proved that their procedure leads to an unconditional convergence to the solution and that the 
truncation error is of order 0(St 2 + SM 2 ) where 6t is the discrete time-step and 6M the grid-step for the field. 

We limited ourselves to the potential Va (<p) = ■^f'-P 2 + 3t5<^ 4 with A = 10 (arbitrary units) and g = 1.2 being kept 
fixed while varying r in order to make some contact with previous numerical studies (2fjj |. This is a well acceptable 
choice of initial condition provided one does not want to describe nonuniversal characteristics of, e.g., the Ising model. 
The model was studied with the sharp, Litim and smooth cut-off regulators introduced in Section [II B[ Some studies 
with the simplified smooth cut-off described by Eq. (|29| with v — 1/2 were also performed to compute x(M) in the 
vicinity of ±Mo. Since no numerical simulations are available for this continuous version of the (p 4 model we tried to 
extract from our data similarities as well as differences between the various versions of the LPA. 

Typically we retained M max = 8, after checking that the results were not influenced by the boundaries conditions. 
The time-step was typically 5t = 2.10~ 4 and the field-step ranging from 5M = 2.10 -3 to SM = 2.10~ 6 according to 
the problem at hand. In all cases the integration of the EDP was stopped at some k m - ln between ~ 10 -6 and ~ 10 -10 
after checking that the solution did not evolve any more. Typically about N t ~ 10 5 time steps were necessary to 
reach convergence. For the smooth cut-off it was moreover necessary to fit the function ui(M) defined as the inverse 
of function C (cf Eq. (|10[) ). its derivatives and inverse yielding unavoidable systematic errors and a reduction of the 
precision on the data. 

1. Thermodynamic potentials below r c 

We display in Fig. [T] the evolution of threshold function Lk (M) when the scale-k of the RG flow decreases from 
A = 10 to fc m in = 10~ 4 . The curves have been obtained with Litim's regulator; other approximations give similar 
curves, except the Sharp cut-off since Lk can become negative in this case. The universal behavior (|2"61 of the LPA 
is exemplified in Fig. [5] in the appropriate reduced variables. 

The behavior of the potential Uk{M) which can be extracted from Lk(M) is displayed in Fig.[3J Note that the second 
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FIG. 10: Left : magnetization versus r c — r for Litim (solid line), sharp (dashed line), smooth cut-off (dashed-dotted line) 
regulators. Right : inverse susceptibility x" 1 = U k=0 (M = 0) versus r — r c , same symbols. In both cases <p 4 model with A = 10 
with g — 1.2 

derivative the true Gibbs potential U k {M) of the k-system remains positive as we shown in sec HID I (cf. Eq. (|28jl). 
For the considered state and at the lowest k = k m - m = 10~ 4 where we stopped the integration of the RG flow one 
finds M = 0.9356 and the inverse susceptibility I7' fe ' m . n (M) jumps from 0.3142 at M = M + SM to 0.4512 10~ 10 at 
M = Mo — SM. Lowering the value of fc m i n doe not change the value of the spontaneous magnetization Mo (for a 
given grid step) while U k . (M) may be given any arbitrary small value for \M\ < Mo- 

The approach to convexity is also well illustrated by the anti-clockwise rotation of the plateau of U k (M), i.e. the 
precursor of the magnetic field, towards an horizontal segment as k — > in agreement with Eq. (|2"?| . Finally one notes 
that all the Uk{M) behaves as parabola inside the coexistence curve, according to Eq. (|26l) . as soon as k is sufficiently 
small. 



2. Continuity /discontinuity of the susceptibility 

Here we check numerically the conclusions of Sec HID] concerning the continuity or discontinuity of the susceptibility 
at M = Mo in the various versions of the LPA. Fig 2] displays the inverse susceptibility U k=0 (M) in the vicinity of Mo 
for the sharp, Litim's, and exponential smooth cut-off in d = 3 dimensions. A relatively small field-step SM = 2.10~ 6 
is required to be convinced of the (dis) continuity of the curves. Although such numerical checks cannot be retained as 
mathematical proofs stricto sensu, the theoretical predictions, i. e. continuity of U k=0 (M) at M = Mo for the sharp 
cut-off and discontinuity for the Litim and smooth exponential regulator are quite well supported by the numerical 
outcomes. 

Figs [U [U [7] and [8] display the curves U k=0 (M) obtained numerically in the case of the simplified smooth cut-off 
threshold function, i.e. C(lu) — K (uj — u; ) -1 / 2 for dimensions of space d = 3.25, d — 3, d — 2.75, and d = 2.5 
respectively. Recall that the behavior of any theory involving a smooth regulator is expected to be of that type. 
In all cases A = 10, g = 1.2 and M max — 8. while the parameter r of the ip 4 potential has been adjusted for a 
spontaneous magnetization Mo ~ 0.05. For each dimension "d" we considered several values of field-step SM ranging 
from SM = 2.10~ 3 to SM = 2.10~ 6 corresponding resp. to a sampling of N M = 4.10 3 to N M = 4.10 6 points in the 
interval (0, Mo) in order to emphasize the numerical difficulty to check this (dis)continuity. As expected (see Sec III Dl ) 
a spectacular change of behavior at d = 3. From continuous at low dimensions (cf the analysis of Sec III D Ij) U k=0 (M) 
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TABLE II: Influence of the field-step SM on the critical parameter r c in all the variations of the LPA considered in this work 
: Sharp, Litim and exponential smooth cut-off. r c was determined by dicothomy as discussed in Sec IIII B 31 Reported data 
concern the <^ 4 potential defined at A = 10, g = 1.2 and r = r c . The RG flow equations have been interrupted at fc m i n = 10 -10 . 
In all cases the time-step is St = 2.10 -4 . We report 10 stable digits, numerical uncertainties only affect the 11th. 



SM 


r c (Litim) 


r c (Sharp) 


r c (Smooth) 


2.10" 3 


-0.1969317598 


-0.2971859571 


-0.4120948764 


2.10~ 4 


-0.1969322088 


-0.2971863642 


-0.4120953249 


2.10" 5 


-0.1969322133 


-0.2971863684 


-0.4120953292 



becomes discontinuous for d > 3. Generally at most a single point seems to survive half the way the two discontinuity 
points whatever the tiny the value of SM considered. In particular note that, exactly at d = 3, Uj!_ (M) "looks" 
indisputably discontinuous at M = ±Mo as anticipated by the analysis of sect HID 21 

3. Determination of the critical point and the critical exponents 

It follows from Eq. (|27|) that the quantity yt — U'^(M — 0)/(fc 2 |wo|) is well suited to discriminate states at a 
temperature above r c from those at a temperature below r c . Indeed, as k — > 0, y^ diverges to +00 at r > r c ( since 
the inverse compressibility U^(M — 0) is finite for r > r c ), while it tends to —1 for subcritical states, ( cf Eq. (|27|) ). A 
dicothomy procedure then yields the precise determination of the critical temperature as exemplifies by FigO which 
displays results only for the ultra-sharp cut-off since in the other variations of the LPA (Litim's or smooth cut-off) 
the bunch of curves is roughly similar and brings no new information. The results for r c are reported in table HU For 
a given approximation they depend slightly on the field-step SM but as soon as SM is as small as 2.10~ 4 we noted 
no effect on r c , up to its 12th digit, by decreasing its magnitude. Note that we can ascertain ~ 8 exact digits on 
the value of r c from the data obtained with SM < 2.10~ 4 . Quite noticeable is the dispersion of the values of r c 

TABLE III: Critical exponents of the magnetization and the susceptibility, resp. /3 and 7, obtained by the numerical experiments 
of Sec IIII B~3l The number(s) in brackets denote the linear regression error on the last digit(s). \r — r c \ was varied from 10 -7 
to 10 -4 . The critical exponent a was obtained by using Rushbrook's equality a + 2/3 + 7 = 2. 





Litim 


Sharp 


Smooth 


p 


0.3327(8) 


0.3532(9) 


0.3352(8) 


7 


1.2768(14) 


1.3292(26) 


1.2781(14) 


a 


0.0578(32) 


-0.0356(44) 


0.0515(32) 



corresponding to different cut-off functions which casts some doubts to the validity of the LPA to predict quantitative 
non universal results. This is an important detail that the non-universal parameters appear to depend strongly on 
the initial value considered at k = A. If one changes something (a cutoff function for example), then the nonuniversal 
characteristics of the system changes. Moreover the MF approximation injected as an initial condition for the flow 
is a too crude approximation. From this respect our results are only but superficially at odds with those obtained 
recently by Dupuis and Machado for the Ising and lattice <p A models where an excellent agreement (within a few 
percents) between the LPA prediction for r c and the MC data was found [29] . However in their work all these authors 
adopt a slightly modified version of the NPRG-LPA in which the effective action at k = A is not taken as the mean 
field result but as the exact one (at this scale) and where integrations over the continuous momenta q is replaced by 
a summation over the vectors q of the first Brillouin zone. Similarly in the domain of the theory of liquids, the works 
reported by Reatto et al. attest good agreement between the LPA and the MC data [1314151 ; in this case the exact 
physics of a reference system (the hard spheres fluid) is injected in the theory. 

However when either the spontaneous magnetization Mo(r) (for r < r c ) or the inverse susceptibility U"(M — 0) (for 
r > r c ) are reported on a graph as functions of r — r c rather than versus r all the curves obtained by means of different 
regulators collapse on a single, approximatively universal one, at least at large scale. This striking observation is 
exemplified in Fig. 1101 In first approximation the effect of the cut-off seems to be a simple shift on r c . If we abandon 
the Sirius point of view and zoom in the vicinity of r ~ r c the various routes yield in fact different behaviors since 
the critical exponents differ slightly. A series of numerical experiments which are resumed in Fig. [TT] allows a rough 
determination of the critical exponents of the magnetization -/?- and the susceptibility -7-. As apparent in Table III 
the values obtained for these two exponents in the case of the Litim and Smooth cut-off regulators are in good 
agreement but they differ quite significantly from those obtained for the sharp cut-off. By passing we note that 
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log 10 ( I r - r c I) 

FIG. 11: Bottom : logarithm of the magnetization versus log 10 (jr c — r\) for Litim (circles), sharp (squares), exponential 
smooth cut-off (diamonds) regulators. The lines are linear regressions of the data. Top : logarithm of the inverse susceptibility 
X _1 = U k=0 (M — 0) versus log 10 (|r c — r|), same symbols. In both cases ip 4 model IN d = 3 with A = 10, g = 1.2, and 
MMax = 8. 



Rushbrook's equality yields a negative a (exponent of the specific heat) in the case of the sharp cut-of, a well-known 
flaw of this approximation. A more stringent discussion will be given in next section where the exponents will be 
calculated with a high precision. 



IV. INTEGRATION OF THE ADIMENSIONNED RG EQUATIONS 



A. Adimensionned RG flow equations 

In order to study fixed point solutions and the spectrum of critical exponents of the LPA we need to work with 
dimensionless potentials and fields. Since [U] — k d and [M] = k~ 1+d / 2 we define the dimensionless magnetization 
M = fc 1_d / 2 M and the dimensionless potential U(M) = k~ d U(M), so that the RG flow Eq. © takes now the form 

d t U = dU + (l-^)MU'-2v d C(U"), (55) 

in which the ' designates a derivative with respect to x. Further simplification can be obtained by introducing finally 
reduced variables. In the general case one defines x = M/^/2v2 and u — U/(2v c i). Litim case will be our exception 
with the choice x = M / yj Avd/d, u = dU /(4iv) and the redefinition valid henceforth C(u)) = 1/(1 +ui). With these 
notations we have now in any case : 



d t u = du + (l-~)xu'-£(u") 



(56) 
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Similarly one introduces the reduced threshold potential l(x,t) = Lk(M) — C(u"(x,t)) which obeys the following 
PDE 

w'(l) d t l = 2w(i) + (1 - ~) xu'{l) /' - I" , (57) 

where we recall that u = C^ 1 is the inverse of function cu C(uj). Eq. (I57[) was obtained by differentiating twice 
Eq. ([56]) with respect to the field variable x. 

B. Fixed points and exponents 

Fixed point solutions u*(x) and l*(x) of Eqs. ([5oT) and (|57| resp. are peculiar solutions of the ordinary differential 
equations (ODE) 

= du* + (l- ^)xu*' -£(u*") , (58a) 

= 2w(P) + (l- ^)xJ{l*)l*' -I*" , (58b) 

respectively. Note that if u* is a solution of Eq. (I58ap then Z*(x) — C o u* (x) is a solution of Eq. (|58bl) ; similarly 
from a solution I* of (|58b|) one builds u* (x) = uj o l*(x) to get a solution of (|58a|) . Because of the assumed Z2 
symmetry we will consider only even solutions u*(x) — u*(—x) of Eq. (I58a[) . Then, since l*(x) = C o u* (x), l*(x) is 
also an even solution of Eq. (|58h>[) . 

Once an even fixed point has been found (analytically or numerically) one investigates the behavior of the solutions 
of Eqs. (|56j) and (|57j) in the vicinity of this fixed point. As usual, we write u(x, t) = u*(x) + h\(x) exp(Ai) and consider 
h\{x) as small. Linearizing Eq. (|56p with respect to h\ leads to the eigenvalue equation : 

= (d — A) h\ + (f — ~)xh' x -C'(u*")hl . (59a) 
Similarly writing l(x,t) = l*(x) + 9\(x) exp(At) and linearizing the full RG flow equation with respect to g yields. 

= g x J (I*) (2 - A) + (1 - ^) x [J {I*) g x + u" (I*) I*' gx} - g'x (59b) 

A first glance, since 

l(x,t) = C(u") - C(u*") +c'(u*")hl(x)exp(Xt) 

the two spectra are identical and the eigenfunctions related by 

gx{x) =£ ou*"(x)hl(x). (60) 

Strictly speaking this is true only if h x (x) is not identically equal to 0. Thus if it turns out that h x (x) is either a 
constant or a linear function of x then eigenvalue A does not belong to the spectrum of the RG operator acting on I. 
As well known, odd and even eigenvectors form two mutually orthogonal linear subsets and will be both considered 
in the sequel. 

C. Trivial solutions 

Such trivial solutions exist whatever the type of fixed point; they are more easily detected on the eigenvalue 
problem (|59a[) . i.e. that attached to the linearization of the RG about u*(x). Let us rewrite Eq. (|58a[) for u* as well 
as the equation for its derivative f* — u* ; one has 

= du* + (1- -)xu*' -C{u*"), (61a) 

= + + _£'(„*") (61b) 

Comparing these equations to the eigenvalue problem (|59al) one readily sees that 
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(i) The constant h\(x) = ho is a trivial even eigenfunction of (|59a|) with the eigenvalue A = d 

(ii) h\ — f*(x) is also an eigenfunction of (|59a[) associated to the eigenvalue A = — 1 + d/2. Since u*(x) should be 
even then f*(x) is an odd eigenvector; according to Wegner it is associated with the shift operator }28| . 

(iii) h\ — x with A = 1 + d/2 is another trivial odd eigenvector associated with the magnetic field. 

From the remark of the previous section it follows that only the trivial odd eigenvalue A = —1 + d/2 survives in the 
eigenvalue problem (|59a[> . is. that attached to the linearization of the RG flow about l*(x). According to (|60p the 
corresponding eigenvector is the odd function gx=-i+d/2(x) = £ ° /* (x) f* (x). 



D. Gaussian fixed point 




We show now that, provided that d > 2, the LPA admits a Gaussian fixed point, with the usual spectrum of 
exponents, irrespective to the type of cut-off. The fixed-point equation for u*, i.e. Eq. (I58a[) . admits u* (x) = 
as a special solution. By integration it gives u* (x) = (Z2 symmetry) and u*(x) = uq with duo — £(0); where 
£(0) = 0, 1 for the sharp and Litim's cut-off respectively and, numerically, £(0) — 5.9973827 for our smooth cut-off. 
Obviously l*(x) — Iq is the related special solution of Eq. (|58b>|) provided Iq — duo- 

The linearized eigenvalue problem for u(x,t) reads 

(d - A) h(x) + (l-^)xh!(x)-£' (0) h" (x)=0, (62) 

where £'(0) < (cf £'(0) = -1 for Sharp and Litim's cut-off, while, numerically £'(0) = -1.37960752 for the 
smooth-cut-off). The change of variables 

y = fix (63a) 

(2-^(0) 2 "=^ («■» 

allows to rewrite Eq. (|62"T) as Hcrmitc equation : 

H"(y)-2yH'(y)+2nH(y) = 0, (64) 

If we request the potential to be bounded by polynomials, n must be restricted to an integer. Then Eq. (|64[) becomes 
the equation defining Hermite Polynomials H n (y). The parity of H n (y) being that of n, the spectrum of even 
eigenvalues is given by 

X p = d - p {d - 2) , (65) 

with n = 2p and p = 0, 1, 2 . . .. The trivial relevant operator Hq with Ao = d is indeed present in the spectrum as 
discussed in sec (|IV C|) . Whatever the type of cut-off we recover the well-known result : in d = 4 we find one relevant 
operator (H2 with Ai = 2), a marginal operator (H4 with A2 = 0) and the first irrelevant operator with A2 = —2 (Hq) 
while in d — 3 there are two relevant operators (H2 and with Ai = 2 and A2 = 1, respectively) and a marginal 
operator (Hq with A3 =0), while the first irrelevant operator is Hg, with A4 = — 1. Some remarks are in order. 

• An analysis similar to that of rcf. Q reveals that the marginal operators become irrelevant at the quadratic 
order. 



The trivial eigenvalue A p= o = d disappears from the spectrum of Eq. (|59b|) . The eigenvectors of Eq. (|59bj) are 
given by g\ p (x) = (4/(2 - d)) H^y) as follows from Eq. ([60]). 



• The odd spectrum is given by A p = d— (p+ 1/2) (tf— 2), p an integer. It includes the trivial solution Ao = l + d/2. 
The other trivial odd eigenvalue A = — l + d/2 is absent accidently from the spectrum, since it should correspond 
to a zero eigenvector (h\ = Jq = 0). 
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FIG. 12: Singular field x c as a function of the initial value 1(0) of Eq. (|58b|) in the case of the Smooth cut-off. 

E. Non Gaussian fixed point 

In this section we focus on the non Gaussian fixed point in d > 2. Recently, the LPA fixed point equation for 
u*(x) have been solved with a very high numerical precision for the sharp and Litim's cut-off (22l.[23|. Here we report 
numerical solutions for l*(x) only in dimension d = 3 ; the three cut-offs considered in this paper were examined and 
compared. Eq (I58b [) can be solved by the shooting method with boundary conditions imposed either at the origin 
x = or at x = oo[22|, (23[, hence the names given to the two methods considered below. 

1. ab origine 



TABLE IV: Results from the shooting method "ab origine" . 





Litim 


Sharp 


Smooth 


r(o) 


0.122859820243702 


0.61903040294652 


7.355923051 


u*"(0) 


-0.186064249470314 


-.46153372011621 


-.7995141985 




20.644305503116 


94.128646935418 


22.95208767 



The LPA fixed point equation (|58bj) is a second order ODE for the function l*(x). Here we solve it numerically 
in d = 3 by providing two initial conditions at x = 0. The first one, I* (x) = ensures the parity x) = l*(x) 
required by Z2 symmetry. It is now well-known that an arbitrary value of ^*(0) yields a solution singular at some 
finite value x c of the field. At x c we have l*(x c ) = for the Smooth or Litim's cut-off and l*(x c ) = —00 for the sharp 
cut-off. These singular values for l*(x c ) correspond to u* (x c ) = +00. Actually, the general solution of (|58b[) involves 
a moving singularity of the form 

l(x) ~ KWd - 2 (x c - x) 1/2 x]J 2 Smooth, Litim (66a) 

l(x) ~ In ((x c ~ x) x c (d/2 - 1)) Sharp (66b) 

where, in the case of Litim's cut-off K\ — 1 (the behavior near the singularity x c is driven by the asymptotics at 
infinity of the function C(u ), cf Eq. (|16l0 . Figure [T2l displays the variation of the singular point x c with the initial 
condition 1(0) in the case of the Smooth cut-off (similar curves are obtained for the Sharp and Litim's cut-off). Two 
peaks where x c diverges to 00 can be noticed. The one on the left corresponds the Gaussian fixed point where l(x) 
is a constant with 1(0) = Ig = £(0) = 5.9973827 (see section llVD|) . The one one the right corresponds to the 
Wilson-Fisher fixed point with 1(0) — 7.355923051 that we are looking at.. 

Our requirement is that the physical solution must be non singular on the entire range x € (0, 00) so we must push 
x c to infinity by adjusting the value of 1(0) by a dichotomy process in the vicinity of the right peak of figure [12] 33]. 
Of course high precision ODE solvers are required for this kind of study. 




400 




FIG. 13: Fixed point solution l^\x) for the smooth cut-off. Red line : shooting ab origine, black line : shooting ad originem. 



We solved equation (|58b[) as well as all the ODE of this paper with the DOPRI853 algorithm of Hairer et al. 
[3lj which is an explicit Runge-Kutta integrator of order "8" and order "5" embedded, with adaptive step-size. We 
imposed relative and absolute errors of 10. ~ 18 and 10. _25 respectively; the code was written in FORTRAN90 in 
quadruple precision. Even with this high technology it is impossible to obtain very large values of x c by tuning 1(0). 
Our results are summarized in Table HVl Our results for u* (0) deduced of our result for l*(0) agree within 14 figures 
with those of ref.[22|, [23| in the case of Litim's and sharp cut-off regulators. The data reported in Table HVl for the 
smooth cut-off are much less precise due to the use of fits for computing function C{uS). The complicated behavior of 
function l*[x) is exemplified in fig. (|T3")) (red curve); We displayed only the result for the smooth cut-off, Litim's case 
is similar while for the sharp cut-off l*(x) tends to — oo at x c instead of "0". 

This ab origine method is rather deceptive but however usefull to check the data obtained by the ad originem 
shooting method that we discuss now. 



2. ad originem 

The existence of a moving singularity at x c which seems impossible to "push" at infinity suggests to solve Eq. ()58b|) 
with initial conditions at infinity, at least a: max large, towards x = 0. The analysis of the asymptotic solutions of 
(I58bl) for x — > oo shows the existence of power law solutions. This second family of solutions with regular scaling 
properties must obviously be preferred to the singular solutions discussed in section (jIVE lj) . One can show that, 
asymptotically, for x — > oo one has, for the smooth cut-off regulator 

t( x )= bl x-^ +^bfx- 2 ?< + 



Ci 



+ (§ + §) bfx-^+0(x-^), 



(67a) 



with /3i — 4/(d — 2) and bi an arbitrary coefficient. Recall that CI, C2, C3 enters the asymptotic behavior (flTf of 
u>(l*) as Z* — > and are given in tableQ] Litim's case can be obtained from (I67a[) by the substitution C% = 1, C2 = —1 
and C3 = 0. In the case of the sharp cut-off one has 

I* (x) = - In bi - Pi In x - -x-P' - 

h 



bid(d- 2) 



x 



- ft - 2 + 0(a;- ft - 4 ) . (67b) 
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FIG. 14: Histogram p(x) of the order parameter of the (d=3) Ising model at the critical point. Black : MC data [34|, red 
LPA (Litim), green : LPA (Sharp), blue : LPA (Smooth) 



In all cases the asymptotics depend on a single parameter bi which fixes the two initial conditions l*(x) and I* (x) 
x = i max . bi is determined in such a way that I* (0) = 0. This is the shooting method ad originem. 
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TABLE V: Results for the coefficient 6j*. 





Litim 


Sharp 


Smooth 


&? 


33.3250777220334 


.091029082436564 


392.879344670467 



The code DOPRI853 detects a stiff problem at large "x" so it is impossible to choose a very large value of x max . 
Actually we retained x max = 15, 22 and 20 for Litim's, the sharp and smooth cut-off regulators respectively. However 
taking into account the full asymptotic expansions ()67a|) or (I67b[) one recovers exactly, i.e. with all the significant 
figures reported in table ITVl . the values of l*(0) obtained by the shooting ab origine. The values of coefficient b* at 
the fixed point obtained by a dichotomy are given in table [V] Stiff integrators could be used instead of DOPRI853. 

Figure [T31 displays our results for l*(x) (smooth cut-off case) obtained by the two shooting methods. At small "x" 
both curves coincide quasi exactly up to some x; at x two branches of the solution separate one that comes from the 
origin, the other, with scaling properties at infinity which comes from infinity. This hysteresis phenomenom cannot 
be discarded numerically. 

From I* one deduce u* = uj(1*) and thus the fixed point u* by integration (up to an additional constant). In 
private and informal discussions the opinion circulates that the histogram of the order parameter p(x) of the 3d Ising 
model at its critical point should coincide with exp(— u*(x)) (up to normalization constants). We are not aware of any 
rigorous proof of this assertion but tentatively took it seriously. This histogram has been obtained by Monte Carlo 
simulations [34[ and is compared in figurdLH with the theoretical predictions of the LPA within the three versions 
considered in this paper. Note that all the histograms has been normalized in such way that the two first even 
moments are equal to unity, i.e. J dxp(x) = 1 and J dxp(x)x 2 = 1. To inject some quantitative elements in the 
discussion of these curves we note that the kurtosis K — J dxp(x)x 4 is K = 1.60399 experimentally, while in the 
LPA one finds K = 1.86162, 1.70294 and 1.51128 for Litim's, the sharp and smooth regulators respectively. 

The behavior of p(x) at large deviations "x" can also be obtained in the framework of the LPA by computing the 
asymptotics of the fixed point solution u*(x). Assuming a power law behavior as x — > +oo, Eq. (|58aj) is used to obtain 



u*(x) 



b* x^ + ^ 



K2 
b* 2 



(d-2) 2 



{d-2f 



2f 



(2d) 2 (d + 2) 4 (d + 4) 



-2(ft,-2) + 



(68a) 



with /3 U = pi + 2 = 2d/(d — 2) and 6* a coefficient which enters at each order in the asymptotics u*(x) (it's value is 
such that u* (0) =0). Of course &* is related to the coefficient 6 ; * which governs the asymptotics of l*(x), one finds 
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TABLE VI: Results from the critical exponents of the LPA. 





Litim 


bnarp 


bmootn 


M 




1.4oU4124ol f f 412 


l.oo f 0o 


V 


0.649561773880648 


0.689459056162135 


0.650594 


m 


0.655745939193339 


(j.oyozoyoozzozoDi 


0.054104 




3.180006512059168 


2.8o842oo582417oo 


O 1 71 OO 

3.17183 


CJ 3 


C HI 07/I77H1 

5. y 122,30012747701 


5.184192105o59884 


5.55112 


W4 




/ .oyo f 92o804o0411 


8. (69 I 2 


CJ5 


1 1 7nono7ctrooopotr7 

11.79oUo765ooo6o57 


10.057908901)049430 


11.7609 


oJe, 


14.8yoODol /ODOOZyo 


12.000 ( 22422o8901o 


1 A O A TO 

14.84 (6 


\ 


n /i 000000000000000 
u.4yyyyyyyyyyyyyyy 


n /i 00000000000000 
u.4yyyyyyyyyyyyyy 


U. 000004 




1.8867038380914204 


1.691338925641807 


1.88197 


ui2 


4.5243907336707728 


3.998514715824934 


4.51219 




7.3376506433543136 


6.382503789820088 


7.31630 


U)4 


10.2839007240259583 


8.821709390384049 


10.2522 


^5 


13.3361699643459432 


11.302996690411253 


13.2933 



that b u = Ki(d — 2) 2 /(b\2d(d + 2)). Recall finally that K±, K2 enters the asymptotic behavior ([TrJf of / = C(uS) as 
u) — > +00 and are given in table HI 

Litim's case can be obtained from (|68ap by the substitution K\ = 1 and K2 = — 1 while, in the case of the sharp 
cut-off, one has 

+ (2 + ln(&«0»(/3u - !))) I d , ( 68b ) 
where, once again /3 U = 2d/(d - 2) but 6,* = b* u 2d{d + 2)/(d - 2) 2 . 



3. critical exponents 

We turn now to the eigenvalue equation (|59bj) . Again this is a second order ODE the solution of which is charac- 
terized a priori by two integration constants. Actually one of these is irrelevant and corresponds to the arbitrariness 
of the normalization of an eigenfunction. Since the fixed point solution l*(x) is an even function of x equation (|59bp 
is invariant under a parity change and the spectrum separates into even and odd eigenvalues. The second integration 
constant is thus fixed by the choice g (0) = (even) or g(0) — (odd). The shooting ad originem method seems 
mandatory and one proceeds as in section IIVE21 

For the smooth cut-off regulator the asymptotic behavior of the eigenfunctions as x — > +00 is found to be 

g(x)~ S(x ai +2^ b*x a '-< 3 > + 

S^gf + g)^- 2 * + ...). (69a) 

where S is an arbitrary constant and on = (A + 2)/(l — d/2). Imposing this form for g{x) at some large x max one 
tunes A to obtain either g (0) = or g(0) = The case of the sharp cut-off is special : 

g(x)~S( X »*-?^ + ...), (69b) 

where oi x = A/(l - d/2). 

Our results for the even and odd spectra are reported in table IVll In the even case there are no trivial eigenvalues, 
as shown in Sec. IIV CI and the first eigenvalue Ai is related to the critical exponent v = 1/Ai of the correlation 
length of the Ising model. The first negative eigenvalue A2 is minus the Ising-likc first correction-to-scaling exponent 
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uj\ — — A2 and so-on. In the odd case, as discussed in Sec. IIV CI there are in general two positive trivial eigenvalues 
Ai and A2, among which only the second one A2 = d/2 — 2 = 0.5 survives in the spectrum. The first non-trivial 
eigenvalue is negative and defines the subcritical exponent 6*5 = uj\ = — A3 and so on. 

The numerical data for A2 = 0.5 serves as a stringent test for the precision of the numerical procedure; while (at 
least) 15 significative digits are obtained for the sharp and Litim's cut-off no more than 6 digits can be ascertain in 
the case of the smooth cut-off. It originates in the various fitting procedures devised to evaluate the function to(l) 
and its derivatives. In the case of Litim's and the sharp cut-off our results are in perfect agreement with those of 
ref. [IH 23] . Overall good agreement between the Litim and smooth cut-off spectra should be stressed. 

Since 77 = in the LPA, one can compute all critical exponents from the critical exponent of the correlation length 
v by the scaling relations. One has for d — 3 : a — 2 — 3^, f3 = v/2, and 7 = 2v. The values are reported in table IVlH 
The comparison with the data of Table IIII1 obtained by solving the dimensionned PDE flow equation, is deceptive 
since there is no good agreement even by taking into account the error bars. The values obtained for (3 ( 7 ) by 
the numerical experiments of section IIIII are systematically larger (smaller) than that obtained in this section. The 
explanation of this discrepancy relies probably in systematic errors due to the size of the field and time steps in the 
numerical resolution of the PDE. 



TABLE VII: Critical exponents in the LPA 





Litim 


Sharp 


Smooth 


p 


0.324780886940324 


0.3447295280810675 


0.325297 


1 


1.299123547761296 


1.37891811232427 


1.301188 


a 


0.051314678358056 


-0.068377168486405 


0.048218 



V. CONCLUSION 

In this paper we attempted to make an exhaustive study of the properties of the ip 4 model in the ordered phase in 
the framework of the NPRG within the LPA approximation. 

We shown that the approach to the convexity is independent of the cut-off, but that fine details are strongly 
affected by the choice of the regulator, notably the analytical behavior of the inverse magnetic susceptibility x~ 1 (M) 
at M — ±Mq. We proved that the singularity of the threshold function C(w) about its largest pole (or essential 
singularity) uq governs the behavior of x _1 (M) at M = ±Mo; if C(uj) ~ (to — coo)~ 1 ' as u — > ljq+ then the inverse 
compressibility is discontinuous at M = ±Mo (as required on physical grounds) only if the inequality 4 — d — 2v < 
holds. In particular Litim's and any smooth cut-off yield a discontinuity of x _1 (Af) in dimension d = 3, while the 
sharp cut-off incorrectly predicts a continuous behavior and thus a merging of the spinodal and binodal curves. 

We have confirmed these subtle properties of the solution below r c by extensive numerical experiments with the 
help of a new algorithm which solves the RG flow equations for the threshold functions rather than for the potential 
or one of its derivatives. The main advantage of the method is to replace the numerical resolution of a highly non 
linear PDE which exhibits numerical instabilities in the ordered phase by that of a quasi-linear parabolic PDE with 
good convergence properties. 

The " standard" version of the LPA retained in this work does not allow to compute the critical parameter r c of the 
model which depends strongly upon the choice of cut-off. The main reason is that the choice of the MF approximation 
as an initial condition for the RG flow is a too crude approximation. Modified versions of the NPRG [l^-LTM Hsl 
remedy to this flaw and yield a quite precise estimate for r c . However we noticed that, choosing as a new variable 
r ~ r c instead of r, the thermodynamics of the LPA (spontaneous magnetization, magnetic susceptibility) for the if 4 
potential is remarkably independent of the cut-off, except very close to the critical point. The latter merely shift r c to 
incorrect values. The numerical solution of the dimensionned RG flow equation does not yield a very precise estimate 
of the critical exponents either, probably because of small numerical errors in the resolution of the PDE. 

In order to compute precisely the critical exponents one must solve the linearized RG about the Fisher fixed point. 
It can also be done for the threshold functions instead of the potential or its derivatives without any noticeable 
numerical differences but with the advantage of getting rid of some trivial solutions corresponding to redundant 
operators. The solution of the resulting fixed point equations and associated eigenvalue problems can be obtained 
with a high numerical precision with the help of a non-stiff solver like DOPRI853 31] for instance. Since the ad 
originem problem is stiff, stiff integrators could be of some help however to reduce the integration step and should be 
tested. 
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It is not clear whether the LPA scenario for the ordered phase survives for more elaborate approximation schemes, 
this could be the subject of further investigations. 
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